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1. INTRODUCTION 

Since modem ultra-spectral remote sensors have thousands of channels, it is difficult to include all of them in a lD-var 
retrieval system. We will describe a physical inversion algorithm, which includes all available channels for the atmospheric 
temperature, moisture, cloud, and surface parameter retrievals. Both the forward model and the inversion algorithm 
compress the channel radiances into super channels. These super channels are obtained by projecting the radiance spectra 
onto a set of pre-calculated eigenvectors. The forward model provides both super channel properties and jacobian in EOF 
space directly. For ultra-spectral sensors such as Infrared Atmospheric Sounding Interferometer (IASI) and the NPOESS 
Airborne Sounder Testbed Interferometer (NAST), a compression ratio of more than 80 can be achieved, leading to a 
significant reduction in computations involved in an inversion process. Results will be shown applying the algorithm to real 
IASI and NAST data. 


2. RAPID RADIATIVE TRANSFER MODEL 

The rapid radiative transfer model used in the lD-var retrieval is a Principal Component-based Radiative Transfer Model 
(PCRTM). The theoretical basis of the PCRTM has been described by Liu et al. [1]. Unlike traditional rapid radiative 
transfer models, which either predict channel radiances or transmittances, the PCRTM predicts the PC scores of the radiance 
spectrum, which contains information for all the channels. The relationship between the PC scores and the predictors, 
monochromatic radiances, is derived from the properties of eigenvectors and instrument line shape (ILS) functions. Each PC 

score (Yj) is calculated via a dot product of the corresponding eigenvector ( U , ) and the channel spectrum ( R chan ) a s shown 
in equation 1 . 

Y i =U i -R chan (1) 

PCRTM is much faster than channel-based RTM since it predicts PC scores (or “Super Channel’’) directly instead of channel 
radiance or transmittance individually. The radiance variation as a function of temperature, FLO, O 3 , CH 4 , N 2 0, CO, surface 
emissivity and observation geometry is captured via monochromatic radiative transfer calculations. The redundant spectral 
information is captured via EOF representation. 

3. METHODOLOGY OF 1D-VAR RETRIEVAL 

The first step in lD-var retrieval is to convert the observed radiance into super channels by projecting the observed radiance 
spectrum onto each individual principal component determined during the PCRTM training process. The current inversion 
algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance and a priori information as 
constraints: 

*„ + i - U = (K r S; l K + W + S-'r'K r S- y , Ky n - yj + K(x„ - x a ) (1) 

where the subscripts n and a represent iteration number and a priori , respectively. Y m is the super channel from observation. 
Y„ is the forward model calculated super channel using the state vector obtained from the nth iteration. S y and S a are error 



covariance matrices associated with y and background state vector x a , respectively, D is a diagonal matrix whose elements 
are determined by the diagonal elements of the k t S> 'k matrix. The Jacobian (or K matrix) is generated by the PCRTM 

model. The dimension of the super-channel used in our retrieval system is on the order of 100, which is about a factor of 84 
less than the original channel radiance spectrum. 

4. RESULTS AND CONCLUSIONS 

We will show results obtained during a Joint Airborne IASI Validation Experiment (JAIVEx), which was held during April 
and May 2007. lD-var retrievals were performed using super channels obtained from the IASI (Infrared Atmospheric 
Sounding Interferometer) and NAST-I (NPOESS Airborne Sounder Testbed Interferometer) instruments. The IASI has 8461 
spectral channels and NAST-I has 8632 channels. It would be very time consuming to perform a lD-var retrieval using all 
8000 channels if a channel-based retrieval approach was taken. However, using the PCRTM retrieval approach, we can 
effective including information from all 8000 channels in the retrieval, while reducing the dimension of data by more than 
one order of magnitude. 
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